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SUMMARY 

The experimentally obtained results of various investigators of 
liquid-metal heat-transfer characteristics were examined and found to be 
not always directly comparable because of differences in experimental 
apparatus or in methods of calculation. The experimental data were 
therefore re-evaluated using as consistent assumptions end methods as 
possible and then compared with each other and with theoretical results. 

The re-evaluated data for both local fully developed and average 
Nusselt numbers in the turbulent flow region were found still to have 
considerable spread, with the bulk of the data being lower than predicted 
by existing analysis. An eqmtion based on empirical grounds w6ich best 
represents most of the fully developed heat-transfer data is 

Nu = 0.625 Pe®’^ 

vhere Nu represents the Nusselt mamber and Pe^ the Peclet nuinher. 

The theoretical prediction of the heat transfer in the entrance region 
was found to give lower values., in most cases, than those foiand in the 
experimental work. 

The theoretical and experimental results for the ratio of local 
Nusselt mmiber to fully developed Nusselt number were integrated to ob- 
tain predictions for the ratio of average Nusselt number to fiilly devel- 
oped Nusselt number for a range of Peclet numbers and length-diameter 
ratios. Most of the e^erimental data faJ_l between 60 to 80 percent of 
the predicted values. 

The experimental evidence was insufficient to serve as a basis for 
any conclusion concerning liquid-metal heat transfer in the laminar or 
transition flow regions. 
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INTRODUCTION 

The use of liquid metals as heat -transfer media is presently of con- 
siderable interest, A number of theoretical and experimental investiga- 
tions to determine the heat-transfer characteristics of liquid met ails 
have been made by various investigators (refs. 1 to 26). In the liter- 
ature^ the results of the experimental investigations often have been 
compared with each other and with the results of theoreticaJ. investiga- 
tions. During the course of investigations of liquid-metal heat-transfer 
characteristics at the NACA Lewis laboratory^ the work of the various 
experimentaJ. investigators was carefully examined. It was found that 
different investigations were not always directly comparable becaiise of 
differences in the experimental apparatus or in the methods of calcula- 
tion. Some of the differences found were: 

(a) Liquid -metal physical properties different from those currently 
accepted were sometimes used. 

(b) At times^ center-line temperatures in and out of the test sec- 
tion were measured rather than "mixing-cup" temperatures. 

(c) Some of the experiments were conducted with unifomn heat input 
to the wall of the test section^ while others more closely approached 
constant wall temperature. 

(d) Some investigators measured the combined heat -transfer coeffi- 
cient in a tube and concentric annulus; different methods were used to 
obtain the individual coefficients. 

(e) Some investigators measured local fully developed heat-transfer 
coefficients; others measured average over-all coefficients. 

(f) The velocity profiles entering the test section varied; some 
approached a fully developed turbulent profile^ while others were more 
nearly uniform. 

(g) Different length-diameter ratios of test section were used. 

The differences in experimental apparatus of items (f) and (g) affect 
only the average heat -transfer coefficient and not the fiilly developed 
coefficient . 

Because of the differences in experimental apparatus and methods of 
calculation listed^ the experimental data of references 1 to 26 were re- 
evaluated using consistent assumptions and methods in order to permit a 
better intercomparison of the experimental results and conparison with 
the resiilts of theoretical investigations. 
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SOBOLS 

a constant 

c specific heat, Btu/(lb)(°F) 

D equivalent or hydraulic diameter, ft 

annulus inner diameter, ft 
Dq annulus outer diameter, ft 

f friction factor 

G weight flow per unit area, lb/(hr)(sq ft) 

Gz Graetz number, PeD/L 

h heat-transfer coefficient, Btu/(hr)(sq ft)(°F) 
k thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 

L length of test section, ft 

m constant, eq. (lO) 

Nu Nusselt number, hD/k 

n constant, eq. (lO) 

Pe Peclet number, RePr, GDc/k 

Pr Prandtl nimnber, c^i/k 

Re Reynolds number, GD/p. 

St Stanton number, h/cG 

t^ fluid center-line temperature, °F 

t fluid buQk temperature, ("builk: temperature" as used in this 

report is synonymous with "mixing-cup tenperature" and "mixed 
mean temperat\ire " ) 

t^ wall tenperature, °F 

X distance along test section, ft 

p fluid bulk viscosity, lb/(hr)(ft) 
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Subscripts and superscripts refer to: 
a average 

f fully developed 

X at station x 

' annulus 


PROCEDURE 

The experimental data of the various references were re-evaluated 
as consistently as possible, plotted as Nusselt number against Peclet 
or Graetz number or against both, and the results compared with theo- 
retical. predictions. These three steps will be discussed in reverse 
order because some of the methods used in re-evailuating the data were 
determined by theoretical considerations. 


Theoretical Investigations of Liquid-Metal Heat Transfer 

The following discussion gives a brief description of some of the 
results of theoretical investigations and is not intended to be cong^lete. 
All the theoretical, investigations discussed consider only the turbulent 
flow region. 

Fully developed heat-transfer coefficients . - Heat -transfer coeffi- 
cients for liquid metals in turb\jlent flow with fully developed velocity 
and temperature profiles have been predicted by a number of investigators 
using somewhat different assumptions. 

(a) Uniform heat input to the wall; rottnd tubes: The most frequently 

analyzed case is that of heat transfer to a round tube with uniform rate 
of heat input along the length of the tube. This case was investigated 
by Martinelli (ref. 2?) using the "momentum transfer analogy." Lyon 
(ref. 6) found a simplified equation which approximated Martinelli 's 
more con 5 )lex relation. This equation, which is recommended by the 
Liquid-Metals Handbook (ref. 28), is 

Nuf = 7.0 + 0.025 Pe^^'® (l) 

Cope (ref. 29) investigated the possibility of assuming that the "modi- 
fied vorticity transfer analogy" applied to the turbulent core of the 
fluid, while the "momentimi transfer analogy" applied to the boundary 
layer and buffer layer. Kennison (ref. 30) assiomed that the heat trans- 
fer is analogous to the transfer of vorticity for turbulent fluid flow 
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in a long straight pipe. Deissler (ref. 3l) modified the "momentum 
transfer analogy" to allow for heat transferred by conduction to or from 
a turbulent particle as it moves radially in the tube. Deissler 's ana- 
lysis is for a Prandtl number of 0.01. 


Some of the resiilts of these various investigations are shown in 
figure 1. The experimental results for fully developed heat transfer 
in a round tube with uniform heat input will be compared with Dyon's 
equation (eq. (l)) inasmuch as this is the equation recommended by the 
Liquid -Metals Handbook and most commonly used in practice. 


(b) Uniform wall ten 5 )eratvirej round tubes: The fully developed 

heat-transfer coefficient in tirrbulent flow in a round tube with a uni- 
form wall temperature has been investigated by Seban and Shimazaki (ref. 
32) using the "momentum transfer analogy"; they give, as an approximate 
relation, the eqixation 


0.8 

= 5.0 + 0.025 Pe^ 


( 2 ) 


This eqixation is plotted in figure 1. The Liquid -Metals Handbook lists 
the equation as 

NUf = 4.8 + 0.025 Pe^*^*® (3) 

and gives the work of Seban and Shimazaki as a reference. The experi- 
mental results for fully developed heat transfer in a round tube with a 
imiform wall temperature will be compared with Seban and Shimazaki 's 
equation (eq. (2)). 


(c) Uniform heat input; annuli: Very little theoretical work has 

been done on the fully developed heat -transfer coefficient in ann\jli. 
For thin anniili (diameter ratio less than or eqxial to 1.4), the Liquid- 
Metals Handbook recommends the use of the theoretical relation proposed 
by Seban (ref. 33) for heat transfer to parallel plates with heat 
through one side only. 

Nu^ = 5.8 + 0.020 (Pep for Dq/d^ 1.4 (4) 

For annuli of diameter ratio greater than 1.4, the Liquid-Metals Hand- 
book lists an equation which approximates the results of Bailey (ref. 
34) and is of the form suggested by Werner, Tldball, and King (ref. 7). 

O r 0 8T 

Nu^ = 0.75 (Do/l>i) ■ 7.0 + 0.025 (Pep ' for D^/D^ >1.4 

( 5 ) 

Equations (4) and (5) are plotted in figure 2. The experimental data 
on heat transfer in anniili will be compared with these eqmtions. 
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Local heat-transfer coefficients in entrance region . - Heat -transfer' 
coefficients in the entrance region have been calculated by several in- 
vestigators for a number of different cases. Poppendiek, Palmer, and 
Harrison (refs. 26, 35, and 36) have analyzed the case of uniform wall 
temperature for various different entering velocity profiles; the ana- 
lysis assxames the eddy diffusivity of heat is negligible when compared 
with the molecular diffusivity and consequently is intended only for low 
Reynolds mmibers. The analysis is independent of Prandtl number. 

Deissler (ref. 37) analyzed the case of uniform heat input at the wall, 
with a f\ally developed velocity profile at the entrance; the n\americal 
calculations were carried out only for a Prandtl number of 0.01. Seban 
and Shimazaki (ref. 38) have made calculations for the case of uniform 
wall ten 5 >erature and fully developed velocity profile at the entrance 
for a Prandtl number of 0.01 and Reynolds numbers of 10^ and 10^. The 
results of the analyses of Poppendiek and Palmer and of Deissler are 
shown in figure 3. 

Average heat-transfer coefficients . - Predictions of average heat- 
transfer coefficients can be made by integrating the predictions for 
local heat-transfer coefficients over the length-diameter ratio of the > 
tube in question. Heat transfer in the entrance region, however, has 
been analyzed for only relatively specialized cases . Therefore, the ex- 
perimental results for average heat-transfer coefficients will first be 
compared with equations (l) and (2), even thou^ equations (l) and (2) 
are derived for fiilly developed heat-transfer coefficients. Later in 
the report, a comparison will be made with the average heat-transfer 
coefficients on the basis of the analytical evidence. 

Teniperature distribution . - The fiilly developed teniperature distri- 
bution due to heat transfer to a liquid metal in turbulent flow in a 
round tube has been predicted on theoretical grounds by several investi- 
gators. The predictions of Martinelli (ref. 27) are shown in figure 4 
for a Prandtl number of 0.022. Martinelli, using his owa predicted val- 
ues for the ten 5 >erature distribution, calculated the ratio of the tem- 
perature differences (t^ - V)/(V - t^,) as a function of Reynolds and 

Prandtl numbers. Martinelli 's results are shown in figure 5 for Prandtl 
numbers pertinent to liquid metals . Martinelli also calculated values 
of (t^ - tjj^)/(t^ - t^) for fully developed flow between flat plates with 
heat flow through both walls with imiform heat flux. These resialts are 
shown in figure 6. 


Methods of Calculation 

The heat -transfer parameters were evaluated using the same method 
of calculation for each individual reference as was used by the authors 
of that particular reference, with the following exceptions: 
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(a) All physical properties of liquid metails were taien from the 
second edition of the Liquid-Metals Handbook (ref. 28). These properties 
are shown in figures 7 to 10. 


(b) When an investigator measured the combined liquid -metal heat- 
transfer coefficient in a tube and concentric annulus, the individual 
heat-transfer coefficients were obtained by assuming that the ratio of 
the Nusselt number in the tube to the Nusselt number in the annvilus is 
determined by eqiiations (l), (4), and (s). 


Nu 

Nu' 


7.0 + 0.025 Pe^'^ 
5.8 + 0.020 Pe'0-8 


for ^q/^± ^ 1-4: 


( 6 ) 


Nu 7.0 + 0.025 Pe^-Q 

0.75(Do/Di)°'^ (7.0 + 0.025 Pe'°'® 

Because of the lack of theoretical work on average heat -transfer 
coefficients^ particLilarly in annuli^ the sajne ratios which have heen 
assumed for the fully developed Nusselt numbers will he assumed for the 
average Nusselt numbers . 

It is important to note thslt in most of those tests in which the 
combined coefficient in a tube and concentric annulus was measured^ the 
Reynolds number in the annulus was smaller than the Reynolds number in 
the tube. Quite often the flow in the annulus was in the trsinsition 
flow region^ while the flow ±n the tube was in the turbulent flow region. 
Inasmuch as there are no predictions for liquid-metal heat transfer in 
the transition region, eqiiations (6) and (?) will be used to separate 
the tube and annulus heat-transfer coefficients even when the flow in 
the annulus is in the transition region. This procedure is open to 
question, and the inteip)retation of the data calculated by this proce- 
dure may be inaccurate. 

(c) In those tests in which the center-line temperature of the fluid 
was measured instead of the bulk temperature, the tenqperature difference 
between the wall and the bulk fluid will be calculated from Martinelli^s 
relation for (t^ - ** "^c^ 5 and 6). Martinelli*s predic- 

tion of (t^ - ■bjj^)/(t^ “ "^c^ flat plates with heat flowing through 

both sides will be used for annuli inasmuch as no other predictions 
covering as broad a range of Reynolds and Prandtl numbers axe available. 



RE-EVAimTION OF EXPERIMENTAL DATA 

The experimental investigations of references 1 to 26 will first be 
discussed individually and then coir 5 )eLred with each other and with theory. 
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The experimental work of the various investigators will be discussed- 
in a chronological order determined by the publication date of the orig- 
inal manuscript. 

Styrikovich and Semenovker . - Styrikovich and Semenovker (ref. l) 
investigated heat transfer to mercury as part of their investigation of 
the mercxiry-steam binary power cycle. Their test sections were a series 
of five tubes, each about 106 inches in length; the diameters were 0.63, 
0.87, 1.58, 1.67, and 1.97 inches. The tubes were heated by external, 
electric heaters. Thermocouples were placed 17.2 inches apart on the 
outside surface of each tube. The b\uLk fluid ten?)erature in the test 
section was calculated by adding to the inlet temperature the temperature 
rise corresponding to the heat input. The velocity profile of the mer- 
cury entering the test section was essentially fully developed. The 
method of heating the mercury approximated uniform heat input to the 
wall. The heat-transfer coefficient was calculated for only the central 
portion of the tube. The coefficients presented are essentially the 
fully developed heat-transfer coefficients . 

The physical properties used in evaluating the heat-transfer coeffi- 
cients are not listed, but the Prandtl number is tabulated over a range 
of temperature from 32° to 1112° F. These Prandtl numbers are lower 
than the values in reference 28, which lists values of Prandtl number up 
to 600° F. The values of Prandtl number of Styrikovich and Semenovker 
and of reference 28 ane shown in figure 11. Since the specific heat and 
viscosity in the temperature range used are essentially the same in ref- 
erence 28 as those reported in the International Critical Tables (l929 
edition), the inaccuracies in Prandtl numiber may be assumed due to in- 
correct values of thermal conductivity. It appears that Styrikovich and 
Semenovker used the thermal conductivity data of Gelhoff and Weumeier, 
which have been found to be high (ref. 12). It was deemed advisable to 
recalculate the data of Styrikovich and Semenovker using the values of 
thermal conductivity from reference 28. The precise temperature level 
of the various data points is not reported, but the average temperature 
level is given as about 932° F. At this temperature, Styrikovich and 
Semenovker list a Prandtl mmiber of 0.CXD56. Reference 28 presents 
Prandtl nimiber data up to 600° F which when extrapolated to 932° F give 
a Prandtl numiber between 0.006 and 0.007. The data points were re- 
evaluated using a Prandtl number of 0.0065 at 932° F. This increased 
the Nusselt and Peclet numibers of the data by about 16 percent. The re- 
evaduated data of Styrikovich and Semenovker are shown in figure 12; 
shown for comparison is equation (l). 

Om i 1 and , Musser , and Page . - Gilliland, Musser, and Page (refs. 2 
and 3) measured both heating and cooling coefficients for mercury. The 
heating test section was 0.319 inch in inside diameter by 14 inches in 
length; heat was added by dropwise condensation of steam on the outside 
of the test section. The cooling test section was 0.319 inch in inside 
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diameter by 51 inches in length; it was cooled by water flowing on the 
outside in a direction opposite to that of the inside flow. Both test 
sections were made of nickel. The mercury and water bulk temperatures 
entering and leaving the test sections and the stream temperature and 
pressure entering the test section were measured. The velocity profile 
of the mercury was ftilly developed at the entrance to both test sections 
The methods of heating and cooling the merciuy were such that the heat- 
ing tests approximated a constant wall temperature;, while the cooling 
tests were somewhere between a constant wall temperature and a consteuat 
heat input. The heat-transfer coefficient measured was an over-all 
average coefficient. . 

Inasmuch as no wall temperataores were measured, it was necessary to 
separate the mercairy heat-transfer coefficients from those of the steam' 
and water. This was done by the Wilson plot method (see refs. 2 smd 3). 

(a) Heating: Tests were run with water in the place of mercury and 

the Wilson plot method was used to determine the combined resistance of 
the steam film and the wall. The range of water flows covered was svif- 
ficiently small and the scatter of the points sufficiently great that 
values of the combined resistance could be chosen ranging from 40 per- 
cent greater to 15 percent smaller than the value selected. An increase 
of 40 percent in steam and wall resistance, however, would increase the 
mercury coefficient only about 8 percent. An attempt was made to use 
the resixlts of the mercury nms to confirm the steam end wall resistance 
but in this case the range of mercury flows and data scatter permit 
selecting a value of resistance ranging from 200 percent greater to 50 
percent lower than the value chosen. The slope of the Wilson plot for 
the runs with water can be compared with the slope predicted by the 
standard empirical, relation for heat transfer to water (ref. 39, p. 168) 

Wu = 0.023 Re®’® Pr®’^ (8) 

The slope predicted by equation (8) turns out to be considerably higher 
than the slope best representing the experimental data. 


(b) Cooling: At a given merciiry flow rate, the water flow rate was 

varied and the combined mercury film and tube wall resistance determined 
by means of a Wilson plot. The range and scatter of the data are such 
that the resistance of the mercury and the wall covild be chosen 20 per- 
cent lower or 15 percent higher than the value actually chosen. The 
corresponding variation in mercury coefficient would be somewhat greater 
Alternatively, cooling coefficients for mercury were calculated by eval- 
\mting the coefficients for water in an annulus using the following 
equation (ref. 39, p. 202): 


St Pr 


2/3 _ 


0.020 


(Do/Di) 


0.53 


Re 


0.2 


( 9 ) 
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The resulting mercury coefficients were approximately 40 percent lower 
than those derived hy the Wilson plot method. The physical properties 
used by Gilliland^ Musser, and Page are about the same as those of ref- 
erence 29. In view of the possible inaccuracies in the method of eval- 
uating the data^ the reported res\ilts of Gilliland, Musser, and Page 
may not be very accurate. Their data are shown in figure 13 without 
change; shown for comparison are equations (l) and (2). The lower val- 
ues for cooling coefficient may be due to the longer length-diameter 
ratio of the cooling section. 

Elser. - Elser (ref. 4) measured cooling heat-transfer coefficients 
for mercury. Three different test sections were, used: The test-section 

inner diameters were 0.317, 0.308, and 0.260 inch; the 0.317-inch- 
■diameter test section was made of mild steel; the other test sections 
were made of stainless steel. The test sections were all over 38 Inches 
long, but measurements were made between two stations 10.2 and 38.3 
inches from the entrance. The mercury was cooled by water flowing in a 
concentric annulus in a direction opposite to the flow of mercury. Two 
thermocouples imbedded in the wall measured the wall temperature at the 
two stations. Two other thermocouples immersed in the stream measured 
a temperature close to the fluid center-line ten^ierature . The velocity 
profile of the mercury at the first station was fully developed. The 
cooling water flow rate was such that a uniform heat input to the wall 
was approximated. The fully developed heat-transfer coefficient was 
measured. 

The only mercury property listed by Elser is Prandtl n\jmber. These 
values are in agreement with the values of reference 28, and the other 
mercury properties will be ass\amed to be conrect. The basic data are 
not presented by Elser. He presents a plot of Stanton number against 
Reynolds momber showing his data points. The values of Stanton number 
have been corrected by Elser to a common Prandtl number by approximating 
the data with a curve of the form 

Nu^ = a Re^- Pr^ (lO) 

He gives no values of n so that it is impossible to return to the 
basic points . 

Elser measured mercury flow by measuring the mercury pressiire drop 
and assuming the following formula for friction factor: 

4f = Re < 80,000 (ll) 

2 log (ReV^f) - 0.8 Re > 80,000 (l2) 

-^4f 

Equation (ll) is due to Blasius; eqmtlon (l2), to Karman. 
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Elser's heat-transfer coefficients are based on the difference 
between wall and fluid center-line ten 5 )eratures . He is not certain of 
the location (depth) of his wall thermocouples ajid states that the dif- 
ference between a midwall and a wall surface location results in shifts 
of heat-transfer coefficients of 4, 7, and 18 percent, respectively, for 
the three tubes of 0.317, 0.308,. and 0.260 inch diameters.. In Elser's 
data, the wall thermocouple is eissumed to be at the wall midpoint. Mar- 
tinelli 's predictions for- the ratio of the tengjeratiire differences 

- ^m^A^w - ^c) fig- 5) were used to change the heat -transfer 

coefficients of Elser so that they would be based on the difference be- 
tween wall and fluid bulk temperatures . This Increased the Nusselt nxam- 
ber about 40 to 60 percent. In this re-evaluation, the wall themocoup- 
les were assumed to be located at the wall midpoint. If the thermo- 
couples were assumed at the wall siirface, the Niisselt numbers wo\ild be 
somewhat increased. Figure 14 shows the re-evalmted data of Elser; 
shown for comparison is equation (l) . 

Bailey, Cope, and Watson . - Bailey, Cope, and Watson (ref. 5) meas- 
ured cooling coefficients for mercury. The test section was a mild steel 
tube of 0.437 inch inner diameter. The central 18 inches of the tube was 
s^irro^Inded by a water jacket, with about 6 inches left projecting at each 
end. These ends were enclosed in chambers in such a manner that the in- 
let and outlet mercury passed along the outside of the ends before enter- 
ing and after leaving the test section. Fluid temperatures were measured 
at the inlet and outlet of the test section; wall temperatures were meas- 
ured at fo\ir stations along the length of the water- jacketed section of 
the tube. 

There is considerable question as to just what tenperatirre was meas- 
\ired at the test-section outlet. First, there was no provision made for 
mixing before the exit temperature was measured. Second, inasmuch as 
the mercury was being cooled, the temperature distribution of the merciury 
was such that the temperature near the wall was lower than the bulk tem- 
peratiire. The merciiry was discharged from the test section into a larger 
chamber, turned 180°, and passed over the end of the test section which 
projected from the water jacket. Because of the mixing in the discharge 
and turning processes, the merciiry on the outside of the projecting end 
had a temperatiure profile close to flat. Hence, the mercury on the out- 
side of the projecting end of the test section was at about fluid bulk 
temperature, while the mercury on the inside of the projecting end (close 
to the wall) was at a temperature lower than fluid bulk temperature. 

Heat was therefore transferred from the outside to the inside; this 
tended to increase the measured mercury exit temperature and consequently 
decrease the observed heat-transfer coefficients. The combined effect on 
heat -transfer coefficient of the heat transferred through the projecting 
end and the lack of mixing before the exit temperature measurement is 
very difficiilt to estimate. 
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The velocity profile at the entrance to the water- jacketed section 
of the test section was close to fully developed. The method of cooling 
was such that uniform wall temperature was approximated at the lower mer- 
cury Peclet nvmibers, while uniform heat input to the wall was more nearly ■ 
the case at hi^ merciiry Peclet mmibers. Fully developed heat-transfer 
coefficients were meas\ired. 

The physical properties used by Bailey, Cope, and Watson were some- 
what different from the values of reference 28j the Prandtl numbers were 
about 10 percent high. The data of Bailey, Cope, and Watson were there- 
fore re-evaluated, using the physical properties of reference 28, in two 
ways: first, it was assiomed that the measiired mercirry exit temperature 

was equal to the fluid b-ulk temperatiire; second, it was assimied that the 
measured mercury exit tengserature was equal to the fluid center-line 
ten^jerature and Majrtinelli's predictions of (t^ - bj^)/(t^ " S) 

were used to calculate the fluid bulk tengierature . The results of both 
methods of computation are shown in figure 15; equations (l) and (2) are 
shown for compaxison. Because of the uncertainties described in the 
measurement of merciary exit temperature, it is difficult to say whether 
either set of data in figure 15 is at all correct. 

Ipron . - I^on (ref. 6) used a tube and concentric annvilus to measure 
the combined coefficient resulting from transferring heat from a sodium- 
potass im alloy (52 percent Na, 48 percent K) flowing in the anniilus to 
the same fluid flowing in the tube. The weight flows in the tube and 
ann\ilus are, necessarily, the same. This type of test section is often 
referred to as a "figure eight" and will be so referred to hereinafter. 
Lyon used four different test sections made of nickel and having the 
following dimensions: 


Test section 

A 

B 

C 

D 

Tube inner diameter, in. 

0.432 

0.703 

0.434 

0.434 

Annulus inner diameter, in. 

.500 

.757 

.500 

.500 

Annulus outer diameter, in. 

.715 

.931 

.684 

.684 

Length, in . 

48 

69 

33 

69 


B ulk fluid temperatures were measured at the inlet and outlet of the 
tube and annulus. The velocity profiles of the fluid entering the tube 
and the annulus were approximately flat (uniform velocity). The figure 
eight test section with counter flow gives approximately constant heat 
input to the wall. The heat -transfer coefficients measured were over- 
all average coefficients. 
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Lyon used physical properties which were somewhat different from 
those of reference 28. The specific heat was about 12 percent higher 
and the thermal conductivity was about 6 percent higher. Use of the 
properties of reference 28 decreases both the Nusselt and Peclet numbers 
about 5 percent. I^on assimjed that the resistances of the walls of the 
foiir test sections were approximately constant, neglecting the differ- 
ences in wall thickness . Uyon did not separate the experimental tube 
and annulus coefficients, but rather calciilated a combined predicted 
coefficient using equation (l) for the tube and an equation approxi- 
mating the results of Harrison and Menke (ref. 40 ) 

Nu^ = 4.9 + 0.0175 (Pep°’® (l3) 

for the annulus. The Liquid-Metals Handbook (ref. 28) mentions eqmtion 
(l3), but prefers equations (4) and (s) for heat transfer in an annulus. 

Icon's data were re-evaluated using the physical properties of ref- 
erence 28 and calculating exactly the resistance of the wall. The over- 
all heat-transfer coefficient was divided into a tube coefficient and an 
annulus coefficient assuming that the Nusselt numbers in the tube and 
annulus axe related as in equations (6) and (?), which are taken from, 
equations (l), (4^, and (s). The use of equations (4) and (s) rather 
than equation (l3) for the annulus results in hl^er annulus heat-transfer 
coefficients and lower tube heat-transfer coefficients for the same over- 
all heat-transfer coefficient. The re-evaluated data of loron are shown 
in figures 16 and 17 j eqixatlons (l), (4), and (s) are shown for con^iarison. 

Untermeyer . - The data of Untermeyer were obtained from unclassified 
material in a classified report. Untermeyer measured heating coefficients 
for a lead-bismuth eutectic with and without magneslxim addition. The 
test section was a steel tube 0.25 inch in inner diameter and 18 inches 
in length. The test section was heated by passing electric ciirrent di- 
rectly through it and the fluid it contained. Wall temperatures and 
fluid inlet and outlet ten^peratures were measured. The velocity profile 
at the test section entrance was closer to flat than to fully developed. 

The method of heating most nearly approximated uniform heat input to the 
wall. Local fully developed coefficients were measxired. 

The physical properties used by Untermeyer are different from those 
of reference 28. The thermal conductivity used by Untermeyer was about 
15 percent low and the volumetric specific heat was about 8 percent high. 

It is difficult to determine from the data whether the heat generated 
directly in the fluid has been subtracted from the total heat input. It 
is eilso dlfficiilt to determine whether a mixing chamber was used in the 
measurement of the fluid bulk ten 5 )erature leaving the test section. Fig- 
ure 18 shows the data of Untermeyer re-evaluated using physlcsuL properties 
from reference 28. 
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Wemer, King, and Tldball . - Wemer, King, and Tidball (ref. 7 and 
unclassified data from a classified report) used a figure eight type of 
test section (tube and concentric annulus with same fluid in both) to 
measure heat-transfer coefficients for a sodium-potassium alloy. Cool- 
ing coefficients were measured in the tube and heating coefficients in 
the annulus. Two test sections having the following characteristics 
were used: 


Test section 

A 

B 

Tube inner diameter, in. ^ 

0.68 

0.70 

Annulus inner diameter, in. 

.75 

.75 

Annulus outer diameter, in. 

1.37 

1.37 

Length, in . 

33.8 

33.8 

Material 

304 Stainless 
steel 

Nickel 


The tests in test section A were all run with a 56 percent sodium, 44 
percent potassium alloy. The tests in test section B were run with both 
56 percent sodiim, 44 percent potassium and 23 percent sodium, 77 per- 
cent potassium alloys . Fluid temperatures were measured at the inlet and 
outlet of the tube and of the annulus. In test section A no provision 
was made for mixing the fluid before measuring the outlet temperatures 
of the tube or the annvilxis, except that the fluid turned one right angle 
bend before each thermDco^Jple . The outlet temperatures measured in ex- 
changer A were therefore somewhere between fluid bulk temperatiore and 
fluid center-line temperatxire, probably closer to fluid center-line tem- 
peratiire. In test section B, mixing baffles were used to mix the fluid 
before measviring outlet tenperatures, and the temperatures measured were 
fluid bulk tenperatures . The velocity profiles of the fluid entering 
the tube and the annulus were essentially flat in test section B, and 
between flat and fully developed in test section A. The figure eight 
test section with co^mterflow gives approximately uniform heat input to 
the wall. The heat-transfer coefficients measured were over-all average 
coefficients . 

Wemer, King, and Tidball used physical properties which were about 
the same aa those of reference 28. However, the relation used to divide 
the over-all heat-transfer coefficient in the test section into separate 
coefficients for the tube and annulus is somewhat different from that 
recommended by the Liquid-MeteO-s Handbook (eq. (7)). 
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The experimental data of Werner, King, and Tidball were re-evaluated 
using equation (?) to separate the over-all heat-transfer coefficient in- 
to tube and annulus coefficients. In addition, the predictions of Mar- 
tlnelli for the ratio of the temperature differences (t^. - “ "^c) 

for the tube and the ann\ilus (figs. 5 and 6) were used to make aldowance 
for. the lack of mixing of the fluid before the outlet thermocouples of 
test section A. The re-evaluated data of Werner, King, and Tidball are 
shown in figures 19 to 22; equations (l) and (s) are shown for comparison. 

Sineath . - Sineath (ref. 8) ran heat-transfer tests with mercury in 
rectangular channels. Sineath 's test section was of the figure elg^t 
type except that instead of a tube and concentric anniilus, he had two 
rectangular channels with one common wall. Heat was added to the mercury 
in one channel and removed from the mercviry in the other. The common 
wall of the two channels was 4 inches high by l/4 inch thick and was made 
of mild steel. The channel gap was l/4 inch and the length, 25 inches. 
Fluid temperatures were measvired at the inlet and outlet of the two 
channels . No attempt was made to provide any mixing of the fluid before 
the outlet temperatures were measured except that the abrupt transition 
from a 4- by l/4-inch rectangular channel to the 3/4-inch pipes which 
carried the fluid to and away from the test section probably resulted in 
considerable mixing. The pipe entered the channels at right angles to 
the direction of flow in the channels; there was no smooth transition 
piece between the pipes and the channels. The fluid tenq)eratures meas- 
ured were probably close to the bulk tenperature . However, the abrupt 
change of section at the entrance to the channels probably caused some 
of the heat-transfer surface to be relatively less effective as a resiilt 
of poor local flow distribution. The figure eight test section with 
counterflow approximated xmiform heat input to the wall. The heat- 
transfer coefficients measured were over-all average coefficients. 

Sineath ran four sets of tests. The first three sets were incon- 
clusive because of experimental difficulties with air entrainment and 
with the deposition of mercuroxis oxide on the wall through which heat 
was being transferred. These problems were partially eliminated in the 
fourth set of runs. There was probably no air entrainment during the 
foiirth set of runs; the wall throu^ which heat was being transferred 
was carefully cleaned at the beginning of the runs but was covered with 
a thin layer of scale, at the end. 

Sineath used physical properties similar to those of reference 28. 

The temperatures in the two channels were sufficiently close that the 
heat-transfer coefficients in both channels could be ass^Jmed the same. 

The data of the fourth set of runs of Sineath are shown unchanged 
in figure 23; shown for conparlson is equation (4). The data of Sineath 
are Undoubtedly lower than they should be as a result of the deposit of 
an oxide film on the heat-transfer surface and of the abrupt change of 
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cross section at the entrance to the channels, which makes a portion of - 
the heat-transfer sirrface ineffective. It is difficult, however, to 
estimate the magnitude of these effects . 

English and Barrett . - English and Barrett (refs. 9 and lO) meas- 
ured heating coefficients for mercury. The test sections were of nickel 
and stainless steel 0.051 inch in inner diameter and 0.059 inch in outer 
diameter. The test section length was approximately 1.9 inches. A cop- 
per coating was "bonded to the outside of the test section; the outer 
diameter of the copper was 0.0825 inch. The test section was heated by 
passing electricity directly through it. The inlet and outlet mercury 
bulk temperatures were measured, as was the outside wall temperature 
along the test section; the voltage distribution along the test section 
was also measured. The velocity profile at the test section entrance 
was fully developed. The method of heating most nearly approximated 
uniform heat input to the wall. English and Barrett measured local heat- 
transfer coefficients along the test section and present the local fully 
developed coefficients for all runs. For one run, the local coefficient 
along most of the tube is presented. 

The physical properties used by English and Barrett are the same as 
those of reference 28 except that the viscosity is slightly high at low 
temperat\ires . This- will probably not affect the fully developed heat- 
transfer coefficients, but the entrance region Reynolds numbers should 
be increased 3 to 4 percent. 

The fully developed heat-transfer coefficients of English and Barrett 
are shown unchanged in figure 24;. equation (l) is shown for comparison. 

The entrance heat-transfer coefficients for the one run presented are 
shown in figure 25; Delssler's predicted cxirves for a Prandtl nvunber of 
0.01 and the same Peclet number range are shown for con^jarison. 

Seban. - Seban (ref. U) measured heat-transfer coefficients with 
lead-bismuth eutectic in two different types of test section. One was 
the figure eight type with a tube and concentric anniolus; the other was 
a copper-coated tube heated by external electric heaters . Only combined 
heat-transfer coefficients for the tube and annulus of the figure eight 
test section are presented in reference 11. Not enough basic data (spe- 
cifically fluid temperatiores ) are presented to separate the tube and 
annulus coefficients . Accordingly, only the electrically heated test 
section will be discussed. The test section was 0.652 inch in inner 
diameter by 48 inches in length. The copper coating was for the purpose 
of containing the wall thermocouples in a region of relatively low tem- 
perature gradient and of smoothing out the nonimiformlties of heat input 
of the external electric heaters. The fluid bulk temperatures were meas- . 
ured at the inlet and outlet of the test section, and the wall tempera- 
t\ires were measured at eight stations along the tube. The velocity pro- 
file was close to fully developed at the entrance to the test section. 
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The method of heating approximates very closely imiform heat input to 
the wall. Local heat-transfer coefficients were measured. The local 
fvilly developed coefficients are presented for all the runs; entrance 
coefficients are presented for a few of the runs. 

The physical properties of Sehan are the same as those given in 
reference 28. Sehan had some trouble with fouling, which caused the 
heat-transfer coefficients to decrease with time. Figure 26 shows un- 
changed the fully developed heat-transfer coefficients of Sehan; equa- 
tion (l) is shown for compaxison. Those points taJcen immediately after 
cleaning have higher heat-transfer coefficients than the others. The 
entrance heat-transfer coefficients presented hy Sehan are shown in fig- 
ure 27; predicted curves of Deissler (see fig. 3) for the same range of 
Peclet nimiber and for a Prandtl ntimber of 0.01 are shown for comparison. 

Trefethen . - Trefethen (refs. 12 and 13) used a figure eight type 
of test section to measure heat-transfer coefficients with mercury. Six 
different tubes, described in the following table, were lised in the tube 
and concentric annulus test section: 


Test section 

A 

B 

C 

D 

E 

F 

Tube inner 







diameter^ in. 

0.711 

0.737 

0.585 

0.523 

0.308 

0.429 

Annulus inner 







diameter^ in. 

.749 

.748 

.629 

.627 

.378 

.500 

Annulus outer 







diameter^ in. 

.874 

.874 

.874 

.874 

.874 

.874 

Lengthy in . 

39.3 

39.3 

39.3 

39.3 

39.3 

39.3 

Tube material 

Stainless 

steel 

Stainless 

steel 

Stainless 

steel 

Copper 

Copper 

Copper 


Heating and cooling tests were run in both the tube euid the annu- 
l\as. Trefethen measured the fluid bulk tenq)erature at the inlet and out- 
let of the tube and annuliis. He also measured the wall temperature of 
the outside of the annulus. The velocity profiles were between flat and 
fuilly developed at the entrance to the test section, probably a little 
closer to flat. The counterflow figure eight test section approximated 
a uniform heat input to the wall. Trefethen presents a fully developed 
heat-transfer coefficient for the central section of his tube (from 10.2 
to 29.4 in. from the entrance). In his calculations he assumed that the 
temperature difference between the fluid in the tube and in the annulus 
renalns the same as the terrq)erature difference at the entrance to the 
tube and annulvis, and that the fluid bulk tengjeratiire gradient eilong the 
length of the tube center section is the same an the temperature gradient 
along the annulus outer wall center section. Trefethen separates the 
tube and annulus coefficients in a manner different from that resulting 
from the i^e of eqvtations (6) and (?). 
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The physical properties used by Trefethen are about the same as 
those of reference 28 for the range of ten 5 )erature covered by his exper- 
iments. (The values of thermal conductivity at high ten 5 >erature (ex- 
trapolated by Trefethen to correct the data of Styrikovich Semenov- 
ker) are lower than those of the Liquid-Metals Handbook by about 6 per- 
cent at 212°, 14 percent at 392°, and 20 percent at 662° F . ) 

The data of Trefethen for fully developed heat-transfer coefficients 
were recalculated losing eqxiations (6) and (?) to separate the coeffi- 
cients of the tube and annulus. The re-evaluated data for the tube are 
shown in figure 28j equation (l) is shown for congiarison. Trefethen 
gives enough data to permit the calculation of over-all average heat- 
transfer coefficients for the tube and annulixs. These coefficients also 
were calculated using equations (6) and (?) to separate the individual 
coefficients in the tube and annulus. These data are shown in figures 
29 and 30; eqioations (l), (4), and (s) are shown for conqparison. 

Doody and Younger . - Doody and Younger (refs. 14 and 15) measured 
both heating and cooling coefficients for mercirry with and without sodium 
additions. The test section was a steel tube 0.493 inch in inner diam- 
eter and 61 inches long. The test section was heated or cooled by water 
flowing in a concentric annulus. Both parallel and counterflow inns were 
made. The annulus was 61 inches long, but the annulus entrance and exit 
were each 6 inches from the ends of the test section; the length of the 
test section between the annulus entrance and exit was thus 49 inches. 
Tube wall temperatures were measured at five stations starting l/2 inch 
downstream of the annulus inlet and ending 3 Inches upstream of the 
annulus outlet (inlet and outlet refer to the parallel-flow case); the 

measurements covered 45^ Inches of the test section. The mercury tem- 
perature was measured at the inlet and outlet of the test section. The 
measurement of test-section exit tenq)erature was made without any pre- 
liminary mixing of the fluid, and the exit tenqjerature measured is 
closer to the fluid center-line temperature than to the fluid bulk tem- 
perature. The mercury velocity profile at the test section entrance was 
probably closer to flat than it was to fully developed. The method of 
heating resulted in a wall condition somewhere between uniform heat in- 
put and imlform wall temperature for the counterflow runs. The parallel 
flow runs resulted in a wall condition \diere the rate of heat input 
varied even more rapidly than for the conditioned uniform wall 
tenq)erature . 

The physical properties used by Doody and Younger axe the same as 
those of reference 28 except for the thermal conductivity. The values 
of thermal conductivity used by Doody and Younger are low by about 1 to 
14 percent in the ten^erature range of the investigation. Because of 
the location of the annulvis entrance and exit as described, some effec- 
tive length of test section between 49 inches and 61 inches must be 
selected. Doody and Yoxxnger used a method due to Sherwood and Petrie 
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(ref. 41) and arrived at an effective length of 56 inches. Since the 
vail temperatures at the ends of the test section vere not measured, 

Doody and Yoimger extrapolated the wall temperature meas\irements to 
cover a length of 56 inches, the "effective" length of their test sec- 
tion. Inasmuch as this extrapolation necessarily neglects end effects, 
the data of Doody and Yoimger represent something between over-all aver- 
age and fully developed heat-transfer coefficients. 

The heat balances of Doody and Younger show deviations as great as 
140 percent, with deviations between 40 and 100 percent being quite 
common. Flow was measured by an orifice, and the orifice calibration 
showed variations as great as 50 percent. The end temperature differ- 
ences between the wall and the fluid fovind by the previously mentioned 
extrapolation were very small, varying from about 0.3° to 8° F, with 
values of 2° F or less being extremely common. Small errors in temper- 
ature measTirement can therefore restilt in large errors in log mean tem- 
peratiore difference. 

Doody and Younger attempted to check their experimental apparatus 
by running heat -transfer experiments with butanol. Unfort\mately, most 
of these data were in the transition region. Some of the data were in 
the laminar flow region, and these data were 25 to 75 percent higher 
than the predictions of the Colburn eqviatlon for laminar flow (ref. 39, 
p. 191). 

In view of the diffiCTilties mentioned, the data of Doody and Younger 
may not be very accirrate. The data of Doody and Younger were re-evalmted 
using the physical properties of reference 28 and the predictions of 
Martlnelli for the ratio (t^ - tm)/(t^ - t^) to detenolne the value of 

the temperature differences between the wall temperature and the fluid 
bulk temperature at the test -section exit. The original data on wall 
temperature are not presented in reference 14 or 15; therefore the ex- 
trapolated end temperature will be iised for the wall temperatvire at th6 
test-section inlet and outlet. The re-evaliiated data of Doody and 
Younger are shown in figures 31 and 32; eqmtions (l) and (2) are shown 
for comparison. 

Lubarsky . - Lubarsky (ref. 16) used a figure el£^t type of test 
section to imeastire heating coefficients in a tube and cooling coeffi- 
cients in a concentric anniilus for lead-bismuth eutectic with and with- 
out magnesi\im additions. The test section was 40.2 Inches long, with a 
tube inner diameter of 0.402 inch, an annulus diameter of 0.50 inch, and 
an annulus outer diameter of 0.625 inch. The fluid bulk temperatiires 
at the inlet and outlet of the tube and annulus were measured. The 
entering velocity profile was approximately flat. The figure eight type 
heat exchanger with counterflow approximated tuaiform heat input to the 
wall. Over-all average heat-transfer coefficients were measvired. 
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Lubaxsky used physical properties which were the same as those of 
reference 28. He used equation (6) to separate the heat-transfer coeffi- 
cients in the tube and the annulus. Lubarsky's data are shown unchanged 
in figures 33 and 34j equations (l) and (4) are shown for comparison. 

Johnson, Hartnett, and Clabaugh (lead -bismuth tests) . - Johnson, 
Hartnett, and Clabaugh (refs . 17 and 18) measured heating coefficients 
for lead-bismuth eutectic in an aluminum-coated tube heated externally 
by electric heaters. The test section is very similar to that of Seban 
(ref. 11 ) described above except that an alimiinum coating has been used 
instead of a copper one. The test section inner diameter was 0.652 inch 
and its length, 48 Inches. The fluid bulk temperature was measured at 
the inlet and the outlet of the test section, and the wall ten^jerature 
was measured at ei^t stations along the wall. The inlet velocity pro- 
file was very close to fvilly developed. The method of heating very 
closely approximated uniform heat input to the wall. Local heat -transfer 
coefficients are presented for both the fully developed region and the 
entrance region. 

Physical properties the same as those of reference 28 were used. 

The data on fully developed heat-transfer coefficients are shown unchanged- 
in figure 35; eqioation (l) is shown for comparison. The data on entrance 
heat-transfer coefficients are shown unchanged in figure 36; predicted 
curves of Deissler (fig. 3(b)) for the same x/D and a Prandtl number 
of 0.01 are shown for con 5 )axison; the predictions of Poppendiek and Palmer 
(fig. 3(b)) for lov Reynolds nuinbers are also shown in figure 36. From 
the local data of Johnson, Hartnett, and Clabaugh, it is possible to 
determine an approximate over-a ll average coefficient by plotting and 
Integrating the local coefficients. The results of this procedure are 
shown in figure 37; equation (l) is shown for con^jarlson. 

Isakoff eind Prey . - Isakoff and Drew (refs. 19 and 20) measiired 
heating coefficients with mercury. The test section was a stainless 
steel tube, 0.127 inch in wall thickness, 1.5 Inches in inner diameter, 
by about 223 inches in length, heated externally by electric heaters . 

The fluid bulk tenq)erature was measured at the test section inlet and 
outlet and the outside wall tenq)erature was measured at seven stations 
along the tube. Velocity and tenqjerature profiles in the fluid were 
measured at three stations along the tube ^x/D = 58, 98, and 138) . The 
entrance velocity profile was very close to flat. The method of heating 
approximated very closely unifom heat input to the wall. Local fully 
developed heat -transfer coefficients were measured at the stations of 
x/D eqxial to 98 and 138. The heat-transfer coefficients measiired at 
x/D equal to 58 are still in the entrance region. 

The physical properties used by Isakoff and Drew are the same as 
those of reference 28. The inside wall ten^jerature was calculated in 
two ways: one was to extrapolate the temperature profile in the fluid 

to the wall; the other was to use the measxired outside wall temperature 
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to calculate the temperature drop through the wall. When this calcula- 
tion was made^ the two inside wall ten^ieratures were found to coincide 
for only three of the total of twelve experimental runs. For the other 
nine runs, the inside wall tenrperatxire calculated from the outside wall 
temperatures was higher than the inside wall tengierature as extrapolated 
from the fluid temperature profile. Fluid and wall tenqjeratures for two 
typical runs sure shown in figure 38. 

This discrepancy between the two methods of deter minin g inside wall 
temperature may be due to inaccuracies in the measurement of outside 
wall temperature. The outside wall temperature was measured by eight 
thennocouples at each station. The ten 5 >erature readings of these ther- 
mocouples varied as much as 25 percent from one another in the high flue 
region. The deviation may have been due to the proximity of the thermo- 
couples to the electric heaters. At any rate, the order of magnitude of 
the variation of the thermocouple readings on the outside wall is as 
great as the magnitude of the differences in temperature resulting from 
the two methods of calcxilatlng inside wall tenqperatirre . It is interest- 
ing to note, however, that the inside wall temperature as calculated 
from the outside wall temperature is greater (for nine cases out of 
twelve) than that extrapolated from the fluid temperature profile. This 
effect is that which would be noted if there were some form of inter- 
facial resistance between the fluid and the tube. However, becaiise of 
the circumferential variation of outside wall ten 5 )erature and because 
three of the runs showed no difference in the inside wall tengjerature 
calculated by the two methods, no conclusions can be reached. 

Figure 39 shows the data of Isakoff and Drew for the fully devel- 
oped heat-transfer coefficient (x/D = 138 )j the coefficient is shown for 
both methods of calculating inside wall ten^jeratxirej equation (l) is 
shown for comparison. The entrance data at x/D eqixal to 58 are shown 
in figure 40j inasmuch as the ordinate in figure 40 is a ratio, both 
methods of calculating the inside wall tenqperature give the same results; 
Deissler's predicted curves (fig. 3) are shown for cong)ar Ison . The out- 
side wall tenqjerattire was measured at a siofficlent nuniber of stations 
along the tube (see fig. 38) to permit the estimate of an over-al 1 aver- 
age heat-transfer coefficient. The average over-all coefficient from 
x/D = 6.4 to x/D = 138, based on inside wall temperatures calculated 
from the outside ten 5 >erature, is shown in figure 41. These average 
coefficients are actually lower than the fully developed coefficients 
because the average outside wall teii 5 )erat\ire at the second station 
(x/D = 98) from the test section entrance is higher than might be ex- 
pected from the other measured tenq)eratures . Whether this might be a 
result of the local teng>erature gradients caused by the external elec- 
tric heaters cannot be determined. 

The temperature profiles in the fluid at x/D eqioal to 138 are 
shown in figure 42 and conqpared with the predictions of Martlnelli 
(ref. 27). 
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Stromqulst . - Stromquist (ref. 2l) measxired heating coefficients 
for mercury with and without sodixm additions. The test section was a 
steel tube heated hy passing electricity directly through the tube. 

The following four different test sections were used: 


Test section 

A 

B 

C 

D 

Inner diameter, in. 

0.380 

0.488 

0.753 

0.787 

Outer diameter, in. 

.753 

1.002 

1.500 

1.501 

Length, in. 

47.25 

48.25 

50-25 

50.25 


Fluid bulk temperatiires were measured at the test-section inlet and out- 
let and the outside wall temperatures were measvired at twelve stations 
along the length of the test section. The entrance velocity profile to 
the test section was somewhere between flat and fully developed. The 
method of heating more nearly approximated imiform heat input thsin it 
did uniform wall temperature. Local heat-transfer coefficients, both 
entrance and fully developed, were measxired. 

Stromquist used physical properties which were the same as those 
of reference 28. Figure 43 shows, unchanged, the fully developed heat- 
transfer coefficients of Stromquist; equation (l) is included for com- 
parison. Figure 44 shows, unchanged, the entrance heat-transfer coeffi- 
cient data of Stromquist; Delssler's predicted curves (fig. 3(a)) are 
shown for conparlson. The predictions of Poppendiek and Palmer (fig. 
3(a)) for low Reynolds nunibers are shown in figure 44(a), 

MacDonald and Qulttenton . - MacDonald and Quittenton (refs. 22 and 
23) measured heating coefficients with sodium. The test section con- 
sisted of a monel tube with a copper jacket bonded to the outside. The 
test section was 0.625 inch in inner diameter by 60.05 Inches in length 
and was heated externally by electric heaters . The purpose of the cop- 
per Jacket was the same as the purpose of the copper coating used by 
Seban (ref. ll) and described in a previous section entitled "Seban." 

The fluid bulk temperature was measxired at the inlet and outlet of the 
test section, and the wall temperatures were measured at eleven stations 
along the test section. The entrance velocity profile was close to 
fxxlly developed if the length of piping in the diagram shown in figure 
1 of reference 23 is to scale. The method of heating very closely 
approximated uniform heat input to the walls . Local heat-transfer 
coefficients were measured; the authors present fxxlly developed heat- 
transfer coefficients for a length of the tube from 47.3 inches to 
54.8 inches downstream of the tube entrance. 

MacDonald and Qulttenton xised the same physical properties as ref- 
erence 28. The data for fxxlly developed heat-transfer coefficients are 
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shown^ ^IQchanged^ in figure 45; equation (l) is shown for comparison. 

The data show a great amount of scatter. Consecutive runs at identical 
Peclet nuinbers auid similar temperatiire levels vary as much as 120 per- 
cent in Nusselt numher, with variations of 30 to 60 percent in consec- 
utive runs being common. In view of this scatter^ the entrance coeffi- 
cients and over-all average coefficients have not been calculated, al- 
though the data were sufficient to male these calculations possible. 

Johnson^ Clabau^^ and Hartnett (mercury tests) . - Johnson, Cla- 
baugh, and Hartnett (ref. 24) measured heating coefficients for merciory. 
The test section was almost identicsil to the test section described pre- 
viously under "Johnson, Hartnett, and Clabaugh (lead-bismuth tests)." 

The test section was an aluminum-coated tube 0.652 inch in inner diam- 
eter and 48 inches in length. The fluid bulk temperatiares were measured 
at the inlet and outlet of the test section; the wall temperatures were 
meas\rred at eight stations along the test section. The entrance veloc- 
ity profile was close to fully developed. The method of heating approx- 
imated very closely tmiform heat input to the wall. Local heat-transfer 
coefficients, both fully developed and entrance, are presented. 

The physical properties used by Johnson, Clabaugh, and Hartnett are 
the same as those of reference 28. Figure 46 shows, iinchanged, the fully 
developed heat-transfer coefficients of Johnson, Clabaugh, and Hartnett; 
equation (l) is shown for comparison. The entrance heat-transfer coeffi- 
cients are shown in figure 47. Deissler*s predicted curves (fig. 3(b)) 
are shown for comparison; the predictions of Poppendiek and Palmer (fig. 
3(b)) for low Reynolds number are also shown in figure 47. From the 
experimentally determined entrance and fully developed heat -transfer 
coefficients, it is possible to determine by integration the over-all 
average coefficient. The res\ilting over-all average heat-transfer co- 
efficients are shown in figure 48; equation (l) is shown for coii5>arison. 

Johnson, Hartnett, and Clabaugh (laminar and transition flow) . - 
Johnson, Hartnett, and Clabaugh (ref. 25)~^ve measiared heating coeffi- 
cients for lead-bismuth eutectic and mercury in the 1 ami nar and trans- 
ition flow regions. The test section used was identical to the test 
sections used in the investigations of lead-bismuth and mercury in the 
turbulent flow region by Johnson, Hartnett, and Clabaugh (see the pre- 
ceding sections entitled "Johnson, Hartnett, and Clabaugh (lead-bismuth)" 
and "Johnson, Clabaugh, and Hartnett (mercury)"). The test section was 
an alimiinum-coated tube 0.652 inch in inner diameter by 48 inches in 
length. The fluid b ulk temper atiares were measured at the inlet and out- 
let of the test section and the wall ten 5 )eratuxes were measured at eight 
stations along the test section. The entrance velocity profile was in 
doubt since the flow was mostly in the transition region. The method 
of heating appi*oximated very closely uniform heat input to the wall. 

Local heat-transfer coefficients, both fxilly developed and entrance, 
were measured. 
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The physical properties xised by Johnson; Hartnett; and Clabaugh are • 
the same as those of reference 28. The fully developed heat-transfer 
coefficients are shown in figure 49; equation (l) is shown for comparison. 
The entrance region heat -transfer coefficients are shown in figirre 50. 
There has been no theoretical work on entreuice region heat transfer in 
the transition flow region; therefore no curves can be shown for purposes 
of con 5 )arison. 

Poppendiek emd Harrison . - Poppendiek and Harrison (ref. 26) have 
measured average heating coefficients with mercury in very short test 
sections. The test section was a small hole along the axis of a copper 
disk of 3- inch outer diameter heated on the outside with water. An un- 
heated starting length was used so that the entrance velocity profile 
was very close to fully developed. Three different test sections were 
used: 


Test section 

A 

B 

C 

Inner diameter; in. 

1/8 

l/l6 

1/8 

Length; in . 

1/8 

l/l6 

l/l6 


The fluid bulk temperatures of the mercury were meas\ired at the inlet 
and outlet of the test section. Wall temperatiares were measured at sev- 
eral radial stations in the test section. The method of heating approx- 
imated constant wall ten 5 >erature . 

The properties used by Poppendiek and Harrison are the same as those 
of reference 28. The over-all average heat-transfer coefficients are 
shown in figure 51. Also shown is a predicted curve of Poppendiek and 
Harrison (ref. 26) for average coefficients. They obtained the curve 
by integrating the local coefficients predicted by Poppendiek and Palmer 
for low Peclet numbers (fig. 3). 

So that the data could be compared with the predictions of Deissler’ 
(fig. 3), the local heat-transfer predictions of Deissler were integrated 
to obtain predictions for average Nusselt number for short length-diameter 
ratios. In figure 52; the data of Poppendiek and Harrison for average 
Nusselt number are can 5 )ared with the predictions of Deissler for the aver- 
age Nusselt number for small length-diameter ratios . 


INTERCOMPARISON OF EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 

The experimental results of the various investigators will be com- 
pared with each other and with theoretical predictions. The arrangement 
of the subjects to be considered will be the same as in the section 
"Theoretical. Investigation of Liquid-Metal Heat Transfer" under 
"PROCEDURE." 
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Fiilly developed heat-transfer coefficients . - Fully developed heat- 
transfer coefficients for the case of xxniform heat input to the wall 
were measured by the following investigators from the group of twenty 
investigations reviewed; Styrikovich and Semenovker; Elserj Ihitermeyer; 
English ajid Barrett; Sehan; Trefethen; Johnson, Hartnett and Clahau^ 
(lead -bismuth); Isakoff and Drew; Stromquist; MacDonald and Quittenton; 
Johnson, Clabaugh, and Hajrtnett (mercury); and Johnson, Hartnett, and 
Clabau^ (laminajr and transition flow) . Curves representing mean lines 
through the data of these various investigators are shown in figure 
53(a); when the amount of scatter of a set of data is so great that no 
mean line can be drawn through it, a cross-hatched area is used to rep- 
resent the data; equation (l) is shown in figure 53(a) for purposes of 
con5)arison. The spread of all the data in figure 53(a) is extremely 
great. However, the following data are not considered; 

(a) Data below a Peclet number of 200; These data, being below a 
Reynolds number of 10,000 and therefore in the trajisltion flow region, 
are not intended to be represented by equation (l). 

(b) Data of Elser and of MacDonald and Quittenton; These are less 
reliable because of the very large scatter of the data. 

(c) Data of Untermeyer; Severe corrosion throiighout the dioration 
of the tests caused large changes in the physical dimensions of the test 
section, ‘as well as possibly contaminating the fluid and the heat- 
transfer surface. 

The Nusselt number of remaining data can be compared with the pre- 
dicted values of eqioation (l) as follows; 


Peclet number 

200 

500 

1000 

2000 

5000 

9000 

Range of ratio of 

0.54 

0.55 

0.57 

0.57 

0.55 

0.50 

measured values 

to 

to 

to 

to 

to 

to 

to predicted 
values 

0.69 

0.75 

0.96 

1.04 

1.13 

1.26 


Another method of comparing the data for the fully developed heat 
transfer is to show on a single plot the actual corrected data of al 1 
the investigators. Figure 53(b) is such a plot. If the same data are 
considered valid in this figure as in figure 53(a), a line given by the 
following eqmtlon would best represent most of the data; 

Nu = 0.625 Pe°‘^ (14) 

This equation is purely empirical and does not in any way suggest that 
the theoretical predictions are faulty. However, inasmuch as there is 
a considerable amo\mt of scatter and since most of the data agree fairly 
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well with this line, it would seem preferable for the designer to use 
equation (l4) xmtil further experiment reduces the uncertainty as to the 
precise values of liquid-metal heat-transfer coefficients. 

Local heat-transfer coefficients in entrance region . - Entrance 
region heat-transfer coefficients have been measured by the following 
investigators from the group of twenty investigations reviewed: English 

and Barrett; Seban; Johnson, Hartnett, and Clabaugh (lead-bismuth); 
Isakoff and Drew; Stromquist; Johnson, Clabaugh, and Hartnett (mercury); 
Johnson, Haxtnett, and Clabaugh (laminar and transition flow); and 
Poppendiek and Harrison. The bulk of the data on heat transfer in the 
entrance region is in the reports of Johnson, Haxtnett, and Clabaugh 
(leaxi-bismuth); Stromquist; Johnson, Clabaugh, and Hartnett (mercury); 
and Poppendiek and Harrison; and is presented in figures 36, 44, 47, 
and 52. 

There is considerable scatter in most of the entrance heat-transfer 
data presented. The predictions of Deissler (fig. 3) agree well with 
the data of Stromquist (fig. 44), but fall slightly low when compared 
with the remaining data (figs. 36, 47, and 52). 

English and Barrett, Seban, and Isakoff and Drew present a small 
amount of entrance region heat-transfer data. As may be seen from fig- 
ure 25, the data of English and Barrett are considerably lower than the 
predictions of Deissler. The data of Seban (fig. 27) agree reasonably 
well with the predictions of Deissler. The data of Isakoff and Drew 
(fig. 40) are considerably higher than the predictions of Deissler, 
which may be in part due to the fact that the entrance velocity profile 
of Isakoff and Drew was very nearly flat, while in the analyses of 
Deissler a fully developed velocity profile is assumed at the entrance. 

Average heat-transfer coefficients . - Theoretical predictions of 
the over-all average heat-transfer coefficient can be made from the in- 
formation on local heat-transfer coefficients. The predictions of 
Deissler for the ratio Nu^/Nu^ were integrated mechanically, and the 
values of the ratio of average Nusselt nijmber to fully developed N\isselt 
nvmiber NUg/Nup are shown plotted a.gainst length-diameter ratio x/d 

for various Peclet numbers in figure 54. 

Values of average Nusselt number were determined from the values 
of the ratio NUg/Nvi^ in figure 54 and the values of fixLly developed 
Nusselt number Nuj. of equation (l) . The resiilts are shown in figure 
55, which gives the variation of average Nusselt number with Peclet num- 
ber for several length-diameter ratios. 

The measured average heat transfer coefficients are described. 
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(a) Uniform heat input to the wall; round tubes; Average heat- 
transfer coefficients in round tubes with constant heat input to the 
walls were measured by the following investigators: I^on; Werner, King, 

and Tidball; Trefethen; Lubarsky; Johnson, Hartnett, and Clabaugh (lead- 
bismuth); Isakoff and Drew; and Johnson, Clabaugh, and Hartnett (mercury). 
Curves representing mean lines through the data of these various inves- 
tigators are shown in figure 56. Also shown is the relation for average 
Nusselt number (l/D = 100 ) from figure 55. The Nusselt numbers of the 
data compare with the predicted values for a length-diameter ratio of 
100 as follows (values below a Peclet nxanber of 20O are not considered 
because they fall in the transition flow region): 


Peclet mmiber 

200 

500 

1000 

2000 

5000 

9000 

Range of ratio of 

0.64 

0.75 

0.54 

0.61 

0.75 

0.70 

measured values 

to 

to 

to 

to 

to 

to 

to predicted 
values 

0.88 

1.03 

1.10 

1.21 

0.83 

0.88 


(b) Uniform wall temperature; roimd tubes: Average heat -transfer 

coefficients in roimd tubes with xmiform wall temperatures or with wall 
conditions somewhere between uniform wall ten 5 )erature and uniform heat 
input were measured by Gilliland, Musser, and Page and by Doody and 
Younger. In figure 57 are ciirves representing mean lines through the 
data of Gilliland, Musser, and Page; the data of Doody and Younger sore 
represented by a cross-hatched area because of scatter. Also shown in 
figure 57 aire the relations for average Nusselt number (l/D = lOO) cal- 
culated from equations (l) and (2) and figure 55. The following data 
are not considered: 

(a) Data below a Peclet number of 200: These data are in the 

transition flow region. 

(b) Data of Doody and Younger: The scatter is large. 

The Nusselt number of the remaining data can be ccmpared with the 
predicted values for an L/D of 100 as follows (values of Nusselt 
number halfway between the values of the two theoretical curves in fig. 
57 will be used for corap6a:ison); 


Peclet number 

500 

1000 

Range of ratio of measiired 

0.63 

0.61 

values to predicted 

to 

to 

values 

0.69 

0.72 
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(c) Annuli: Average heat-transfer coefficients in annuli or be- 

tween flat plates with constant heat input to the wall were measured hy 
lyon; Werner, King, and Tidball; Sineath; Trefethen; and Luharsky. Fig- 
ixre 58 shows cinrves representing mean lines through the data of these 
various investigators; also shown are the relations for average Nusselt 
nuniber (l/d = 100 ) calculated from equations (4) and (s) and figure 55. 
The Nusselt numbers of the data compare with the predicted Nusselt num- 
bers (average of the Nusselt numbers of the two theoretical curves of 
fig. 58) as follows (values below a Peclet number of 200 are not con- 
sidered because they fall in the transition flow region): 


Peclet nijmber 

200 

500 

1000 

Range of ratio of measured 

0.69 

0.39 

0.39 

values to predicted 

to 

to 

to 

values 

1.21 

1.37 

0.67 


Temperature distribution . - The only experimental data on temper- 
ature distribution is the work of Isakoff emd Drew. Plots of the tem- 
peratiore distributions measvired by Isakoff and Drew are shown in figure 
42; Martinelli's predicted temperature distributions are shown for com- 
parison. It is possible to use Isakoff and Drew's temperature and ve- 
locity profiles (the measured velocity profiles check quite well with 
the predicted velocity profiles) to calculate the values of the ratio 
(tw ” ” ^c^ shown in figure 59. Martinelli's predictions (fig. 

5) for the ratio (t^ - tjj^)/(t^ “ also shown. The measured val- 

ues are smaller than the predicted values. The predicted values of 
Martinelli for (t^ - tjj^)/(‘t^ ” ^c) we^ vised to calculate the fluid bulk 

temperature in those cases in \diich fluid center-line temperature was 
measured (Elser; Doody and Younger; Werner, King, and Tidball (test 
section A)). If the values of (t^ - tja)/(ty^ " actually lower 

than predicted by Martinelli, the Nusselt numibers of these cases would 
increase . 

Final comparison of heat-transfer data . - The experimental data of 
all the investigators for fully developed and average heat -transfer 
coefficients have been compared with the appropriate prediction and the 
results are shown in table I. 

Figure 60 shows the variation of the ratio of measured to predicted 
Nusselt nvmiber with Peclet number for some of the data of table I. The 
results vdiich are not shown in figure 60 were not included for the fol- 
lowing reasons: 

(a) There is large scatter of data. 

(b) Obvious vincertainties exist as to the accuracy of the data. 
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(c) The measurements are of average heat-transfer coefficients 
which were made conciorrently with the measurement of fully developed 
coefficients; the fxally developed coefficients are shown in figure 60. 

(d) The measurements are of annuliis heat-transfer coefficients 
which were made concurrently with the measxirement of ro\md tube coeffi- 
cients; the rovtnd tube coefficients are shown in figure 60. 

On the basis of the results shown in table I and figure 60, it can 
be seen that most of the measured values of fully developed and average 
Nusselt numbers for turbulent flow (as given by eqs. (l), (2), (4), and 
(s) and fig. 54) fall between 60 to 80 percent of their predicted 
values . 


SUGGESTED EXPERIMENTAL WORK 

It is siiggested that the type of e3q>eriment most likely to reduce 
the uncertainties with respect to liquid-metal heat trsmsfer would be 
one in which velocity and tenperatvLre profiles were measured in the 
fluid, somewhat like the experiment of Isakoff euid Drew. The experiment 
of Iseikoff and Drew could be improved by the use of a thick, high con- 
ductivity, mfetallic coating around the test section simi.lar to the one 
used by Seban or Johnson, Hartnett, and Clabaugh; this woijild probably 
eliminate the loncertainties as the measurements of outside wall 
tenperatxxre . 

The experimental data are insufficient to lead to any conclusion 
concerning liquid-metal heat transfer in the laminar and transition 
flow region. Such data are greatly needed becaiise the small amount of 
data in these flow regions disagree considerably with theoretical 
predictions . 


CONCLUDING REMAEKB 

The review of the experimental investigations of liquid-metal heat 
transfer may be summarized as follows: 

(1) The ejqperimental data of the various investigators were re- 
evalmted using as consistent assumptions and methods as possible, and 
the results were compared with each other and with theory. 

(2) The re-evaluated experimental data for fully developed Niisselt 
number in the tiorbulent flow region were fo\md still to have considerable 
spread, and most of the data are lower than predicted by theory. 
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( 3 ) An equation based on en 5 )irical grounds, which best represents 
most of the fxilly developed heat-transfer data, is 

Nu = 0.625 Pe°’^ 

where Nu and Pe represent Nusselt number and Peclet nimaber, 
respectively. 

( 4 ) The theoretical, predictions of heat transfer in the entrance 
region were found to give lower values, in most cases, than those found 
in the experimental work. 

( 5 ) Integrating the theoretical and experimental results for the 
ratio Nu^/Nu^ gave predictions for the value of the ratio NUg/Nu^ 
over a range of Peclet mmiber and length-diameter ratio. 

(6) The small amount of data on temperatiire distribution disagreed 
with the theoretical predictions, the discrepancy increasing with de- 
creasing Reynolds number. 

( 7 ) The experimental evidence is insufficient to lead to any con- 
clusion about liquid-metal heat transfer in the laminar and transition 
flow region. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 4, 1954 
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TABLE I. - COMPARISON OF HEAT-TRANSFER DATA 


Investigation 

Reference 

Type of heat- 
transfer coefficient 
measured 

Theoretical equation 
used for comparison 

Ratio of measured Nusselt 
number to predicted Nusselt number 
for Peclet niimber cf 

200 

500 

1000 

2000 

5000 

9000 

Styrlkovlch 

and Semenovker 

1 

Round tube; fully devel- 
oped; uniform heat 
Input 

NUf = 7.0 + 0.025 Pej. 



0.80 

0.77 



Gilliland, 

Musser, and Page 
(heating data) 

2, 3 

Round tube; over-all 
average; uniform wall 
temperature 

Muf = 5.0 + 0.025 Pef 
^corrected for L/D = 44 


0.69 

0.72 




Gilliland, Musser, 
and Page (cooling 
data) 

2, 3 

Round tube; over-all 
average; between uni- 
form heat Input and 
uniform wall tempera- 
ture 

NUf = 6.0 + 0.025 Pej. 
(average of eqs. (l) 
and (2)) corrected for 
L/D =160 


0.63 

0.61 




Elser 

4 

Round tube; fully devel- 
oped; uniform heat 
Input 

Nuf * 7.0 + 0.025 PCf 




0.17- 

.40 



Bailey, Cope, and 
Watson 

5 

Round tube; fully devel- 
oped; between uniform 
heat Input and uniform 
wall temperature 

NUf = 6.0 + 0.025 Pef 0-8 
(average of eqs. (l) 
and (2)) 

0.36 

0.46 

0.39 

0.35 



Lyon (tube data) 

6 

Round tube; over -a 11 
average; uniform 
heat Input 

NUf = 7.0 + 0.025 PCf 0-8 
corrected for L/D =110 

0.63 

0.72 

0.80 




Lyon (annulus data) 

6 

Annulus; over-all aver- 
age; uniform heat 
Input 

NUf = 5.8 + 0.020(P6f )0.8 
corrected for L/D = 225 

0.74 

0.90 





Untermeyer (without 
magnesium addi- 
tions ) 


Round tube; fully devel- 
oped; uniform heat 
input 

NUf = 7.0 + 0.025 Pef0*8 


0,16 

0,23 

0,48 



Untermeyer (with 
magnesium addi- 
tions) 


Round tube; fully devel- 
oped; uniform heat 
input 

NUf = 7.0 + 0.025 Pej. 0.8 


0.92 

0.80 




Werner and King 

(heat exchanger A; 
tube data) 


Round tube; over-all 
average; uniform heat 
Input 

NUf = 7.0 + 0.025 PSf 0.8 
corrected for L/D = 49 

0.69 






Werner and King 

(heat exchanger Aj 
annulus data) 


Annulus; over-all 

average; uniform heat 
Input 

Nul = 0.75 (D^j^)®*^ 

[7.0 + 0.025(Pe^)‘^*®] 
Dq/Dj^ = 1.83 corrected 
for L/D = 55 

0.90 

I 

1 

1 

1 

1 

Werner, King, and 
Tldball (heat 
exchanger B; tube 
data) 

7 

Round tube; over-all 
average; uniform heat 
Input 

Nuf = 7.0 + 0.025 Pef 0.8 
corrected for L/D = 49 

■ 



1.14 

■ 

■ 

Werner, King, and 
Tldball (heat 
exchanger B; 
annulus data) 

7 

Annulus , over-all 
average; uniform 
heat Input 

Nu| = 0.75 (Do/Di)°- 5 
[7.0 + 0.025(Pe^)°'^ 
Dq/Dj^ = 1.83; cor- 
rected for L/D = 55 

1.15 

1.29 

1 

1 

1 

■ 

Slneath 

B 

Rectangular ducts; over- 
all average; uniform 
heat Input 

Nuf = 5.8 + 0.020(Pef)0-® 
corrected for L/t) = 50 

■ 

0.41 


■ 

■ 

■ 

English and Barrett 

9, 10 

Round tube; fully devel- 
oped; uniform heat 
Input 

NUf = 7.0 + 0.025 PCf 0.8 , 

0.61 

0.74 

■ 

■ 

■ 


Seban 

.11 

Round tube; fully devel- 
oped; uniform heat 
Input 

NUf = 7.0 + 0.025 Pef 0.8 

■ 

■ 

0.67 

0.58 



Trefethen (fully 
developed tube 
data ) 

12, 13 

Round tube; fully devel- 
oped; uniform heat 
Input 

NUf = 7.0 + 0.025 Pef 0.8 , 

0.68 



0.76 



Trefethen (over- 
all average 
tube data) 

12, 13 

Round tube; over-all 
average; uniform heat 
Input 

NUf = 7.0 + 0.025 Pe^ 0.8 , 
corrected for L/D = 65 

0.74 

0.83 

0.86 

0.87 




^''Corrected for L/b" means that the fully developed Nusselt number found from the equation was multiplied by 
the ratio Nu^/Nu^ from fig. 54. 
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TABLE I. - Concluded. COMPARISON OF HEAT-TRANSFER DATA 


Investigation 


Trefethen (over-all 
average annulus 
data) 


Doody and Younger 
(data with no 
sodium additions) 


Doody and Younger 
(data with sodium 
additions) 


Lubarsky (tube data) 


Lubarsky (annulus 
data) 


Johnson, Hartnett, 
and Clabaugh 
(lead-bismuth; 
fully developed 
data) 


Johnson, Hartnett, 
and Clabaugh 
( lead -bismuth; 
over-all average 
data) 


Type of heat- 
transfer coefficient 
measured 


12, 13 Annulus; over-all 
average; uniform 
heat input 


14, 15 Round tube; cver-all 

average; between uni- 
form heat Input and 
uniform wall temper- 
ature 

14, 15 Round tube; over-all 

average; between uni- 
form heat Input and 
uniform wall temper- 
ature 


16 Round tube; over-all 
average; uniform 
heat input 


Annulus; over-all 
average; uniform 
heat input 


17, 18 Round tube; fully devel 
[ oped; uniform heat 


17, 18 Round tube; over-all 
average; uniform 
heat input 


Theoretical equation 
used for comparison 


Ratio cf measured Nusseit 
number to predicted Nusseit number 
fcr Pec let number of 



Average of Nuf = 5.8 + 
0.020 Pej*®-® and Nuf = 
0.75 (Do/Di) 0-3 I^.O + 
0.025 (Pe^)O-®] , 2, 

corrected for L/D = 20o 

Nuj. = 6.0 + 0.025 Pej. 0.6 
(average of eqs. (l) 
and (2)), 

corrected for L/D = 114 

NUf = 6.0 + 0.025 Pej. 0-8 
(average cf eqs. (l) 
and (2)), corrected for 
L/D =114 


^ = 7,0 + 0.025 P 
corrected fcr L/D 


Nu^ = 5.8 + O.O2O(Pe^)®-0 
corrected fcr L/D » 320 



Isakoff and Drew 

19, 20 

Round tube; fully devel- 

Nu,. = 7.0 + 0.025 Pe 

(fully developed 


oped; uniform heat 


data; inside wall 
temperature calcu- 
lated from fluid 
temperature pro- 


Input 



Isakoff and Drew 
(fully developed 
data; inside wall 
temperature calcu- 
lated from outside 
wall temperature) 


Isakoff and Drew 
(over-all average 
data; inside wall 
temperature calcu- 
lated from outside 
wall temperature) 


MacDonald and 
Qulttenton 


Johnson, Clabaugh, 
and Hartnett 
(mercury; fully 
developed data) 


Johnson, Clabaugh, 
and Hartnett 
(mercury; over- 
all average data) 


Johnson, Hartnett 
and Clabaugh 
(laminar and 
transition flow) 


19, 20 iRound tube; fully devel- INUp = 7.0 + 0. 
I oped; uniform heat j 





Round tube; cver-all 
average; uniform 
heat input 


Nu^ = 7.0 + 0.025 Pe^ 0.6 
corrected for L/D » 138 



21 Round tube; fully devel- Nu- = 7.0 + 0.025 Pe_ ^*8 
oped; uniform heat 
input 


22, 23 Round tube; fully devel- Nu- = 7.0 + 0.025 Pe„ ®*8 
oped; uniform heat 
input 


24 Round tube; fully devel- Nuj. = 7.0 + 0.025 Pej. 0.68 0.70 0.71 0.70 0.65 0.60 

oped; uniform heat 
input 


24 Round tube; over-all Nu = 7.0 + 0.025 Pe 0.8 0.70 0.73 0.75 0.76 0.76 0.69 

average; uniform , 

heat input corrected for L/D = 74 


25 (Round tube; fully devel- INu = 7.0 + 0.025 Pe 0.68 


oped; uniform heat 
input 



































































Peclet number, Pe 


Figure 1. - Theoretical predictions of fully developed Nusselt numbers for heat transfer to 
liquid metals in turbulent flow in round tubes. 
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Figure 2. - Theoretical predictions of fully developed Nusselt numbers for heat transfer to 
liquid metals In turbulent flow In annuli and between flat plates. 
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Figiire 3. - Theoretical predictions of heat transfer to liquid metals in turbulent flow in round tubes. 
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Figure 4. - Theoretical predictions of Martinelli (ref. 27) 
for fully developed temperature profiles for heat transfer 
to liquid metals in round tubes. Prandtl number, 0,022. 
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Figure 6. - Theoretical predictions of Martinelll 
"between flat plates. 
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Figure 7. - Physical properties of mercury (ref. 28). 
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Figure 8. - Physical properties of sodium- potassium alloy (ref. 28). 
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Figure 10. - Physical properties of sodium (ref. 28). 
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Pec let number, Pe 


Figure 12. - Re-evaluated data 
transfer to mercury in round 


of Styrikovich and Semenovker (ref. 
tubes. 


l) for fully developed heat 
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Figure 14. - Re-evaluated data of Elser (ref. 4) for fully developed heat transfer to mercury 
in round tubes. 
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Figure 15. - Re-evaluated data of Bailey, Cope, and Watson (ref. 5) for fully developed heat transfer to mercury in round 
tubes . 
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(a) Test section A (length- diameter ratio, 111). 

Figure 16. - Re-evaluated data of Lyon (ref. 6) for average heat transfer to a sodium- potassium 
alloy in round tubes. 
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Peclet number, Pe 

(a) Test section A 1.43) length- diameter ratio, 223). 

Figure 17. - Re-evaluated data of Lyon (ref. 6) for average heat transfer to a sodiim- 
potassium alloy in annuli. 
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(b) Test section B (D^/d^, 1.23; length- diameter ratio, 397). 

Figure 17. - Continued. Re-evaluated data of Lyon (ref. 6) for average heat transfer 
to a sodium- potassium alloy in annuli. 
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Pec let number, Pe 

(c) Test section C (D^/D^, 1.37; length-diameter ratio, 179). 


Figure 17. - Continued. Re-evaluated data of Lyon (ref. 6) for average heat transfer 
to a sodiiam- potassium alloy in annuli. 
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Figure 19. - Re-evaluated data of Werner, King, and Tidball for average heat transfer to 
sodium- potassium in round tubes. Test section A; length-diameter ratio, 50. 
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Figure 20. - Re-evaluated data of Werner, King, and Tidball for average heat transfer to 
sodium-potassium in annuli. Test section length- diameter ratio, 54; ratio of outer to 

inner diameter, 1.83. 
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Peclet number, Pe 


Figure 21. - Re-evaluated data of Werner, King, and Tidball (^ref. 7) for average heat 
transfer to sodium- potassium in round tubes. Test section B,- length- diameter ratio, 48. 
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Figure 22. - Re-evaluated data of Werner, King, and Tldball (ref. 7) for average heat 

transfer to sodium- potassium in annuli. Test section B; length- diemeter ratio, 54; ratio 
of outer to inner diameter, 1.83. 
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Peclet number, Pe 

Figure 23. - Data of Sineath (ref. 8) for average heat transfer to mercury in rectangular 
ducts. Length-diameter ratio, 50. 




Pec let number, Pe 

Figure 24. - Data of English and Barrett (refs. 9 and 10) for fully developed heat transfer 
to mercury in roiond tubes. 






Figure 25. - Data of English and Barrett (ref. 9) for entrance region heat transfer to mercury in round tubes. Peclet 
number, 800 to 900. 
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Figure 28. - Re-evaluated data of Trefethen (refs. 12 and 13) for fully developed heat 
transfer to mercury in round tubes. 
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Figure 30. - Re-evaluated data of Trefethen (refs. 12 and 13) for average heat transfer to 
mercury In annuli. Length-diameter ratio, 79 to 334. 
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Figure 32. - Re-evaluated data of Doody and Younger (refs. 14 and 15) for heat transfer to 
mercury with small sodium additions in round tubes. 

I 

j 

















Nusselt number, Nu 


I 


78 


NACA TN 3336 



Figure 36. - Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 18) for entrance region 
heat transfer to lead-bismuth eutectic in round tubes. 
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(b) x/D = 13.8. 

Figure 36. - Continued. Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 18) for 
entrance region heat transfer to lead-bismuth eutectic in round tubes. 




Peclet number, Pe 
(c) x/D = 23. 


Figure 36. - Concluded. Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 18) for 
entrajice region heat transfer to lead-bismuth eutectic in round tubes. 
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Figure 37. - Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 18) for average heat 
transfer to lead-bismuth eutectic in round tubes. Length- diameter ratio, 74. 
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Distance along test section, ft 
(a) Run 3. 

Figure 38. - Test-section temperatures from two typical runs of Isakoff and Drew (refs. 19 and 20). 
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Figure 38. - Concluded. Test*section temperatures from two typical runs of Isakoff and Drew (refs. 19 and 20). 
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Figure 40. - Data of Isakoff and Drew (refs. 19 and 20) for entrance region heat transfer to mercury in round tubes (inside 
wall temperatures extrapolated from fluid temperature profile). 
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Figure 41. 
in roimd 


- Data of Isakoff and Drev (refs. 19 and 20) for average heat transfer to mercury 
tubes (inside vail temperatures calculated from outside vail temperatures). 
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Figure 42. - Data of Isakoff and Drew (refs. 19 and 20) for fully 
developed temperature distribution for heat transfer to mercury in 
round tubes. 
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Figure 43. - Data of Stromquist (ref. 2l) for fully developed heat transfer to mercury, vith and without sodium additions, in round tubes. 
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(a) 170 < Pec let number < 1000, 


Figure 44. - Data of Stromquist (ref. 2l) for entrance region heat transfer to mercury, with and without sodium additions, 
in round tubes. 
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(b) 1000 < Pec let number < 2000. 

Figure 44. - Continued. Data of Stromquist (ref. 2l) for entrance region heat transfer to mercury, with and without 
sodium additions, in round tubes. 
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Figure 44. - Continued. Data of Stromquist (ref. 2l) for entrance region heat transfer to mercury, with and without 
sodium additions, in round tubes. 
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□ Pure mercury (nonvetting) 

O Sodium amalgam (vetting) 

O Sodium amalgam (nonvetting) 
Deissler (ref* 37) 
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X 

(d) 3000 < Peclet number < 5000. 

Figure 44. - Continued. Data of Stromquist (ref. 2l) for entrance region heat transfer to mercury, vith and without 
sodium additions, in round tubes. 
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Figure 44. - Concluded. Data of Stromquist (ref. 2l) for entran 
sodium additions, in round tubes. 
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Figure 45. - Data of MacDonald and Quittenton (refs. 22 and 23) for fully developed heat 
transfer to sodium In round tubes. 
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Figure 46. - Data of Johnson, Clabaugh, and Hartnett (ref. 24) for fully developed heat 
transfer to mercury In round tubes. 
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Figure 47. - Data of Johnson, Hartnett, and Clabaugh (ref. 24) for entrance region heat 
transfer to mercury in round tubes. 
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(b) x/D = 13.8. 

Figure 47. - Continued. Data or Johnson, Hartnett, and Clahaugh (ref. 24) for entrance 
region heat transfer to mercury in round tubes. 
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(c) x/D = 23. 

Figure 47. - Concluded. Data of Johnson, Hartnett, and Clabaugh (ref. 24) for entrance 
region heat transfer to mercury In round tubes. 
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Figure 48. - Data of Johnson, Clabaugh, and Hartnett (ref. 24) for average heat transfer 
to mercury in round tubes. Length-diameter ratio, 74. 




Peclet number, Pe 


Figure 49. - Data of Johnson, Hartnett, and Clabaugh (ref. 25) for fully developed heat 
transfer to lead-bismuth eutectic and mercury in round tubes in laminar and transition 
flow regions. 
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Figure 50. - Data of Johnson, Hartnett, and Clabaugh (ref. 25) for entrance region heat 
transfer to lead-bismuth eutectic and mercury in round tubes in laminar and transition 
flow regions. 
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Figure 51. - Data of Poppendlek and Harrison (ref. 26) for average heat transfer to mercury 
in very short round tubes. 
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Figure 52. - Data of Poppendiek and Harrison (ref. 26) for average heat transfer to 
mercury In very short round tubes compared with prediction of Delssler (ref. 37). 



Average Nusselt number, Nu 


104 


MCA TN 3336 



Graetz number, Gz = -j— 

(b) 1000 < Peclet number < 2000. 


Figure 52. - Continued. Data of Poppendiek and Harrison (ref. 26) for average heat transfer 
to mercury In very short ro\md tubes compared with prediction of Deissler (ref. 37). 
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Johnson, Hartnett, and Clabaugh (mercury and 
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(a) Faired curves of experiments. 


Figure 53. - Comparison of measured 
tubes with constant heat input to 


and predicted fully developed -Nusselt numbers in ro\md 
wall . 
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Figure 53. - Concluded. Comparison of measured and predicted fully develo] 
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Figure 54. - Predictions of Deissler (ref. 37) for variation of ratio of average Nusselt number to fully developed Nusselt 
number with length-diameter ratio for various Peclet numbers (Prandtl number^ O.Ol). 
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Figure 55. - Variation of average Nusselt numiber with Peclet nuniber for several length- 
diameter ratios as determined from figure 54 and equation (l). 
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Figure 56. - Comparison of measured and predicted average Nusselt numbers in round tubes 
with uniform heat input to wall. 
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Figure 58. - Comparison of measured and predicted average Nusselt numbers in annuli and 
between flat plates with uniform heat input to wall. 
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Figure 59. - Comparison of measured and predicted values of ratio (t^ - in round tubes vith uniform heat input 

to wall. 
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Figure 60. - Variation of ratio of measured Nusselt number to predicted Nusselt number with Peclet number. 
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